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Abstract. The random walk to be considered takes place in the 5- 
spherical dual of the group U(n + 1), for a fixed finite dimensional irre- 
ducible representation 8 of U(n). The transition matrix comes from the 
three term recursion relation satisfied by a sequence of matrix valued 
orthogonal polynomials built up from the irreducible spherical functions 
of type 5 of SU(rz + f). One of the stochastic models is an urn model 
and the other is a Young diagram model. 



1. Introduction 

Around 1770 D. Bernoulli studied a model for the exchange of heat be- 
tween two bodies. This model can also be seen as a description of the dif- 
fusion of a pair of incompressible gases between two containers. This model 
was independently analyzed by S. Laplace around 1810, see the references 
in [Fj. Another model of similar characteristics was introduced by P. and 
T. Ehrenfest in 1907 in connection with the controversies surrounding the 
work of L. Boltzmann in the kinetic theory of gases dealing with reversibil- 
ity and convergence to equilibrium. Boltzmann had apparentlly deduced 
his H-theorem dictating convergence to equilibrium starting from the time 
reversible equations of Newton. For a nice account of this see [K] . Both of 
these models are instances of discrete time Markov chains with fairly explicit 
tridiagonal one-step transition probability matrices which are obtained by 
considering carefully the underlying stochastic mechanism that connects the 
state of the system at two consecutive values of time. 

The second model features two urns, I and II, that share a total of N 
balls. The state of the system at time n is the number of balls in urn I. 
Each ball has a different label from the set 1, 2, N. At time n a number 
j in the set 1,2, ...,N is chosen with equal probabilities and the ball with 
this label is moved from the urn where it sits to the other urn. This gives 
the state of the system at time n+ 1. Writing down the one-step transition 
probability matrix is now a matter of counting carefully. 

While it had been possible to obtain interesting answers for these two 
models for quite some time, it is only much more recently that some very 
nice connections have been noticed between these models and some basic 
sets of discrete orthogonal polynomials, namely the Krawtchouk and the 
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dual Hahn polynomials. Moreover although there are many ways of arriving 
at these polynomials it is relevant to mention here that they can be realized 
as the "spherical functions" for certain finite bihomogeneous spaces. A very 
good reference for this material is [Sj. We stress the remarkable fact that 
these two models of old vintage and clear physical significance can be solved 
in terms of the simplest of all hyper geometric functions, namely 2F1 and 

3-^2- 

As many readers certainly know many of the classical special functions of 
mathematical physics, such as the Legendre, the Hermite and the Laguerre 
polynomials, could have been obtained for the first time as spherical func- 
tions for certain symmetric spaces. A good basic reference here is [Vj. The 
way that things developed historically is, of course, completely different. 

The interplay between important physical problems and certain tools that 
arise naturally in group representation theory constitutes the theme of this 
paper. The situation described here is the reverse of what has been discussed 
above for the Bernoulli-Laplace and the Ehrenfest models: we will go from 
group representation theory to some concrete models that might be of some 
physical interest. We will start from a matrix that is obtained from group 
representation theory and try to build a model that goes along with it. 
The models constructed here are certainly not the only possible ones. More 
natural ones might be lurking around. 

In a series of papers including [TH HH iGPTl IGPT11 IGPT21 IGPT31 E 
IPTll IPT21 IPT31 IPT4| one considers matrix valued spherical functions as- 
sociated to a pair (G, K) arriving at sequences of matrix valued polyno- 
mials of one real variable satisfying a three term recursion relation whose 
semi-infinite block tridiagonal matrix is stochastic, i.e. the entries are non- 
negative and the sum of the elements in any row is 1. This matrix depends 
on a number of free parameters that have a very definite group theoretical 
meaning. The important point is that the tools developed in the papers just 
mentioned allow one to give explicit expressions, in terms of some definite 
integrals, of all the entries of any power of the original matrix. This means 
that if one could think of a nice Markov chain with this matrix as its one-step 
transition probability matrix one would have an explicit form for the entries 
of the n-step transition probability matrix. Many readers will recognize that 
this is exactly what S. Karlin and J. McGregor, see [KMcG], proposed as 
a way of exploiting orthogonal polynomials and the role they play in the 
spectral analysis of certain finite or semi-infinite tridiagonal matrices. The 
method advocated in [KMcGJ starts with a so called birth-and-death process 
whose one-step tridiagonal transition matrix is easily constructed from the 
given model and one has to look for the corresponding spectral information: 
the eigenfunctions and the spectral measure. Here we travel this road in the 
opposite direction in a more elaborate set-up. 
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The relation between matrix valued orthogonal polynomials, block tridi- 
agonal matrices and Quasi-Birth and Death processes has been first ex- 
ploited independently in [DRSZj [G] as well as in later papers by these au- 
thors. 

We will consider several random walks whose configuration spaces are 
subsets of U(n + l)(k), the so call k-spherical dual of U(n + 1), and whose 
one-step transition matrices come from the stochastic matrix that appears 
in [PT2] and [P], see also [PT4] . The dual of U(n + 1) is the set U(n + 1) 
of all equivalence classes of finite dimensional irreducible representations of 
U(n + 1). These equivalence classes are parametrized by the n + 1-tuples of 
integers m = (mi, . . . , m n+ \) subject to the conditions vn\ > • • • > m n +i. 

If k = (ki, ■ ■ ■ , k n ) S U(n), the k-spherical dual of U(n + 1) is the subset 
U(ra + l)(k) of U(n + 1) of the representations of U(n + 1) whose restriction 
to U(n) contains the representation k. Then it is well known, see [V], that 
U(n + l)(k) corresponds to the set of all m's as above that satisfy the extra 
constraints 

(1) rrii > ki > mj+i, for all i = l,...,n. 

In other words U(n + l)(k) can be visualized as the subset of all points 
m of the integral lattice Z n+1 in the set 

[fcl,Oo) X [k 2 ,ki] X ••• X [kn-Xjkn] x (-oo,k n ]. 

An example is given in the figure below. 



Figure 1. U(n + l)(k), n = 1, k x = 3. 



We can now state more precisely the point of this paper: starting from the 
stochastic matrix M that appears in [PT2] and jP], we describe a random 
mechanism that gives rise to a Markov chain whose state space is the subset 
of U(n + l)(k) of all m 6 U(n + l)(k) such that s m = Sk and k n > 
(s m = mi + • • • + m n+ i, Sk = k\ + • • • + fc n ), and whose one-step transition 
matrix coincides with the one we started from. The construction in |GPT| 
and [PT2] deals with the case of (SU(3), U(2)) but in [PT3] and [F] this was 
extended to the case of (SU(n + l),U(n)). 

One step of the Markov evolution will consist of two substeps taken in 
succesion. In the first substep one of the values of rrii increases by one, 
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subject to the constraints (1). In the second substep one of the new values 
of our mj's decreases by one, again this is subject to the same constraints. 
Thus from the configuration m one could for instance go to m — ej + e,- or 
one could stay put at m. We use the notation ej for the vector with its ith 
component equal to 1 and all the others equal to 0. Any state has a total of 
at most n(n + 1) + 1 positions where it can move in one complete step of our 
process consisting of two simpler steps. Keep in mind that the two succesive 
simpler steps can end up with our random walker in the initial state. We 
will analyze in detail the simpler substeps that constitute one full step of 
our process. This will take up most of the analysis in the next sections. 
We now describe the contents of the paper. 

In Section [2] we collect the necessary material to state and explain a three 
term recursion relation (with matrix coefficients) for a sequence of matrix 
valued orthogonal polynomials, built up from irreducible spherical functions 
of a fixed type associated to the pair (SU(n + l),U(n)). This should help 
the reader make the connection between [PT21 IF] and the present paper. 

In Section [3] we construct a factorization of the stochastic matrix that 
define the three term recursion relation for the sequence of matrix valued 
orthogonal polynomials given in the previous section. This factorization into 
two stochastic matrices leads to the two substeps mentioned above. 

Before starting the analysis of our general urn model in Section [5] for one 
of the substeps, we describe in detail in Section d] an urn model for n = 2. 

The definition of the stochastic matrix M alluded above, as well as its 
factorization make sense for any m 6 U(n + l)(k). 

To each configuration m\ > ni2 > • ■ • > m n > of n integer numbers 
we associate its Young diagram, a combinatorial object which has mi boxes 
in the first row, rri2 boxes in the second row, and so on down to the last 
row which has m n boxes. For example the Young diagram associated to the 
configuration 6>4>4>3is 



Figure 2. 

Young diagrams and their relatives the Young tableaux are very useful in 
representation theory. They provide a convenient way to describe the group 
representations of the symmetric and general linear groups and to study 
their properties. In particular Young diagrams are in one-to-one correspon- 
dence with the irreducible representations of the symmetric group over the 
complex numbers and the irreducible polynomial representations of the gen- 
eral linear groups. They were introduced by Alfred Young in 1900. They 
were then applied to the study of the symmetric group by Georg Frobenius 
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in 1903. Their theory and applications were further developed by many 
mathematicians and there are numerous and interesting applications, be- 
yond representation theory, in combinatorics and algebraic geometry. 

If we consider the subset all m E U(n + l)(k) such that m n +i > it is 
natural to represent such a state of our Markov chain by its Young diagram, 
see Section [H Then in the last two sections we describe a random mechanism 
based on Young diagrams that gives rise to a random walk in the set of all 
Young diagrams of 2n + 1 rows and whose 2j row has kj boxes 1 < j < n, 
and whose transition matrix is Mi, see (I24p . 

2. Spherical functions of (SU(n + 1), U(n)) 

Let G be a locally compact unimodular group and let K be a compact 
subgroup of G. Let K denote the set of all equivalence classes of complex 
finite dimensional irreducible representations of K; for each 5 € K, let £5 
denote the character of 5, d(5) the degree of 5, i.e. the dimension of any 
representation in the 6, and \& = d(8)£,s- We choose the Haar measure dk 
on K normalized by J K dk = 1. We shall denote by V a finite dimensional 
vector space over the field C of complex numbers and by End(V) the space 
of all linear transformations of V into V. 

A spherical function $ on G of type 5 £ K is a continuous function on G 
with values in End(y) such that 

i) ^(e) = /. (i= identity transformation). 

ii) 3>(x)<J>(y) = J K X8{k~ l )$(xky)dk, for all x,y G G. 

If $ : G — > End(y) is a spherical function of type 5 then §{kgk') = 
§(k)$(g)§(k'), for all k,k' G K, g £ G, and k H> 3>(fc) is a representation 
of K such that any irreducible subrepresentation belongs to 5. 

Spherical functions of type 5 arise in a natural way upon considering 
representations of G. If g 1— > U(g) is a continuous representation of G, say 
on a finite dimensional vector space E, then 

P(S) = [ xsik-'W^dk 

JK 

is a projection of E onto P(5)E = E(5). The function : G — > End(£7(5)) 
defined by 

<%)a = P(8)U{g)a, geG,ae E(5) 
is a spherical function of type 5. In fact, if a £ E(5) we have 

$(x)$(y)a = P(8)U(x)P(8)U(y)a= [ xs{k- 1 )P(5)U(x)U(k)U(y)adk 

JK 

= ^xs{k' 1 )^(xky)dk^ a. 

If the representation g i-> C(<7) is irreducible then the associated spherical 
function $ is also irreducible. Conversely, any irreducible spherical function 
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on a compact group G arises in this way from a finite dimensional irreducible 
representation of G. 

The aim of this section is to collect the necessary material to state and 
explain a three term recursion relation for a sequence of matrix valued or- 
thogonal polynomials, built up from irreducible spherical functions of the 
same type associated to the pair (SU(n + 1), S(U(n) x U(l))). 

The irreducible finite dimensional representations of SU(n+l) are restric- 
tion of irreducible representations of U(n + 1), which are parameterized by 
(n + l)-tuples of integers 

m = (mi,m 2 , ■ ■ ■ ,m n+ i) 

such that mi > m 2 > • • • > m n +i. 

Different representations of U(n + 1) can restrict to the same represen- 
tation of G = SU(n + 1). In fact the representations m and p of U(n + 1) 
restrict to the same representation of SU(n + 1) if and only if rrii = pi + j 
for all % = 1, . . . , n + 1 and some j G Z. 

The closed subgroup K = S(U(n) x U(l)) of G is isomorphic to U(n), 
hence its finite dimensional irreducible representations are parameterized 
by the n-tuples of integers 

k = (ki, k 2 , ...,k n ) 

subject to the conditions k\ > k 2 > • • ■ > k n . 

Let k be an irreducible finite dimensional representation of U(n). Then 
k is a subrepresentation of m if and only if the coefficients ki satisfy the 
interlacing property 

mi > ki > mj-fi, for all i = 1, . . . ,n. 

Moreover if k is a subrepresentation of m it appears only once. (See [VKJ). 

The representation space of k is a subspace of the representation space 
V m of m and it is also if-stable. In fact, if A £ U(n), a = (det^) -1 and 
v £ Vk we have 

\0 a) V a \ l) V a \0 l) U ' 

where s m = mi + • • • + m n+ i and Sk = k\ + • • • + k n . This means that the 
representation of K on obtained from m by restriction is parameterized 
by 

(2) (ki + s k - s m , • • • , K + s k - Sm)- 

Let $ m ' k be the spherical function associated to the representation m of 
G and to the subrepresentation k of K. Then ([2]) says that the i^-type of 

$ m < k isk+( Sk -S m )(l,...,l). 

Proposition 2.1. The spherical functions <J> m,k and <& m ,k of the pair 
(G, K) are equivalent if and only if m' = m + j(l, . . . , 1) and k' = k + 
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Proof. The spherical functions $ m,k and $ m ' k ' are equivalent if and only 
if m and m' are equivalent and the -ftT-types of both spherical functions are 
the same, see the discussion in p. 85 of [Tlj . We know that m ~ m' if and 
only if 

m' = m + . . . , 1) for some j G Z. 
Besides, the K types are the same if and only if 

ki + Sk — «m = k[ + Sk' — s m i for alH = 1, . . . , n. 
Therefore k' = k + p(l, . . . , 1), and now it is easy to see that p = j. □ 

The standard representation of U(n + 1) on C n+1 is irreducible and its 
highest weight is (1,0,... ,0). Similarly the representation of U(n+1) on the 
dual of C n+1 is irreducible and its highest weight is (0, ... , 0, —1). Therefore 
we have that 

C n+1 



Vn 



'(1,0,-, o) and ( C T = ^(0,...,0,-l)- 
For any irreducible representation m of U(n + 1) the tensor product V m ® 
C n+1 decomposes as a direct sum of U(n + l)-irreducible representations in 
the following way 



(3) 
and 

(4) 



v m <g> c n+1 ~ y m+ei e v t 



m+e 2 



V, 



m+e n+ i i 



V m ® (C 



n+l\* 



^m— ei © ^111—62 



Vr 



m-e n+ i 5 



where {ei, • • • ,e n+ i} is the cannonical basis of C n+1 , see [VKL 

Remark. The irreducible modules on the right hand side of (J3j) and @ 
whose parameters {m' 1 ,m! 2 , ■ ■ ■ ,m' n+1 ) do not satisfy the conditions > 
m' 2 > • • • > m' n+1 have to be omitted. 

Starting from ([3]) and the following theorem is proved in [P]. 

Theorem 2.2. Let <p and tp be, respectively, the one dimensional spherical 
functions associated to the standard representation of G and its dual. Then 

n+l 

<P(g)^ m 'H9) = J> 4 2 (m,k)$ m + e *> k ( 5 ) 



71+1 



^)$ m - k (5)=2^(m ) k)$ 



m— e;,k 



(g)- 



i=i 



The constants aj(m, k) and 6j(m, k) are given by 



(5) 



aj(m,k) 
6i(m,k) 



mi 



j + i) 



W=i(kj ~ mi-j + i) 



1/2 



mi 



3+i) 



1/2 
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Moreover 

n+l n+1 

(6) J>f(m,k) = ^&?(m,k) = l. 



=1 i=i 




Our Lie group G has the following polar decomposition G = KAK, where 
the abelian subgroup A of G consists of all matrices of the form 



(7) 



(Here I n -\ denotes the identity matrix of size n — 1). Since an irreducible 
spherical function <3? of G of type 5 satisfies <&{kgk') = <&(k)<k(g)<k(k') for 
all k, k' E K and g E G, and <!>(£;) is an irreducible representation of K in 
the class 5, it follows that $ is determined by its restriction to A and its 
K-type. Hence, from now on, we shall consider its restriction to A. 
Let M be the group consisting of all elements of the form 




m= B , B E U(n — 1). 



Thus M is isomorphic to U(n — 1) and its finite dimensional irreducible 
representations are parameterized by the (n — l)-tuples of integers 

t = (£l,t2! • • • > tn-l) 

such that t\ > t2 > • • • > tn-l- 

If a £ A, then <J? m,k (a) commutes with $ m,k (m) for all m £ M. In fact 
we have 

$ m > k (a)$ m ' k (m) = $ m ' k (am) = $ m ' k (ma) = $ m ' k (m)$ m ' k (a). 

The representation of U(n) in Vk C V^, k = (k\, . . . , k n ) restricted to 
U(n — 1) decomposes as the following direct sum 

(8) V k = V t , 

teM 

where the sum is over all the representations t = (ti, . . . ,in-i) £ such 
that the coefficients of t interlace the coefficients of k, that is k{ > ti > /cj+i, 
for all % = 1, . . . , n — 1. Since each Vt C appears only once, by Schur's 
Lemma, it follows that <3? m ' k (a)|y t = <^™' k (a)Id|y t , where (/>™' k (a) G C for 
all a G A. 

By using Proposition 12.11 given a spherical function $ m,k we can assume 
that Sk — s m = 0. In such a case the K-type of <£ m ' k is k, see ([2]). Now it 
is easy to see that if (m, k) is one of such a pair then 

(9) m = m(w,r) = (w + fa, n + k 2 , . . . ,r n -i + k n ,-(w + n-\ hr n _i)), 
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where < w, k n > — (w + r\ + • • • + r n -i) and < rj < ki — ki + \ for 
i = 1, ... n — 1. Thus if we assume w > max{0, — k n } and < ri < ki — ki+\ 
for i = 1, ... n — 1 all the conditions are satisfied. 

We observe that the representations t of M appearing in the right hand 
side of ([8]) are of the form t = r + k', where k' = (&2, • • • , k n ) and r is in the 
following set 

n = {r = (n, . . . ,r n _i) : < r, < fcj - k i+1 }. 
In particular the number of M-modules in the decomposition of Vk is 



N 



n-l 



k 



i+l 



+ !)• 



We will identify <3? m ' k (a) with the column vector (<£ 



m,k / 



plex valued functions $™' k (a) indexed by Q, where <I ) ? 1 ' k (a) 



(a)) re n of N com- 

,m,k / \ 
<Pr+v( a h 
a £ A. 

From now on we fix k G K and take m = m(w, r) as in ([9]) for all 
w > max{0, — k n } and red. Also in the open subset {a(6) G A : < 9 < 
7r/2} of A, we introduce the coordinate t = cos 2 (9) and define on the open 

interval (0, 1) the complex valued function F r!S (w,t) = $™^' r ' ) ' k (a(0)) and 
the corresponding matrix function 

F(w,t) = (i ? r,s('W,t))(r,s)enxn- 

For each w > max{0, — k n } we also define the following matrices of type 
fix!! 



(10) 
where 

(A w ) 



{C u 



A,, 



((A 



B,, 



((B 



w Jr,s J i 



a 2 l+1 (m(w,r))bf (m(w,r) + e n+ i) 
a 2 j+1 (m(w,r))bl(m(w,r) + e j+1 ) 




= ((Cu 

if s 
if s 

otherwise 



r 



a 2 (m(w,r))b 2 n 



■ n+n m(w,r) +ei)) 
al(m(w,r))b 2 +1 (m(w,r) +ei) 




if s 
if s 



r 
r 



(B v 



otherwise 

(m(w , r))b 2 (m(w , r) + ej)) if s = 



l<i<n+l 

(m(u>, r))bl +1 (m(u>, r) + e J+ i 



4+i( m ( w > r ))^+i( m ( w, 1 r )+ e ).+i) 
a| +1 (m(w,r))&f +1 (m(w,r) + e i+ i) 



if s 
if s 
if s 



r + e-i 



r + e,- 







where 1 < i,j < n — 1, and a?(m(io, r)) = 
6f((m + e 3 -,k)) for 1 < z, j < n + 1, see ©. 



otherwise 
af(m,k), 6, 2 (m(w,r) + e,-) 
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Theorem 2.3. For each fixed K-type k = (k\, . . . k n ), for all integers w > 
max{0, — k n } and all < t < 1 we have 

(11) tF(w, t) = A w F(w - 1, t) + B w F(w, t) + C w F(w + 1, t). 

Proof. This result is a consequence of Theorem 12.21 and of the appropriate 
definitions of ^4^, B w , C w given in (llOh . when we take g = a(9). 

We recall that (p(g) and ip{g) are the one dimensional spherical functions 
associated to the G-modules C n+1 and (C n+1 )*, respectively. A direct com- 
putation gives 

4>(a(9)) = (a(9)e n+1 ,e n+1 ) = cos 9. 

and 

tp(a(9)) = (a(0)X n+1 ,X n+1 ) = cos (9. 
Then 0(a(0))V>(a(0)) = cos 2 (#) =t. □ 

If g E G = SU(n + 1) let A(g) denote the nxn left upper corner of g, and 
let A be the dense open subset of all g £ G such that A(g) is nonsingular. As 
in |PT3| in order to determine all irreducible spherical functions of G of type 
k = (k\, . . . , k n ) an auxiliary function : A — > End(Vk) is introduced. It 
is defined by &k(g) = ir(A(g)) where it stands for the unique holomorphic 
representation of GL(n, C) corresponding to the parameter k. It turns out 
that if k n > then <£ k = 3> m ' k where m = (ki, . . . , k n , 0). 

Then instead of looking at a general spherical function $ u, ' r = $ m ( u, ' r )> k 
of type k we look at the function H w,r (g) = & w ' r (g)&k(9) 1 which is well 
defined on A. 

As before we construct the matrix function 

H(w,t) = (# r ,s(?M))(r,s)enxn- 
where H TjS (w,t) = Hs' r (a(9)), t = cos 9 £ (0, 1). 

Let ^(t) = (^r,s(i))(r,s)enxn be the transpose of H(Q,t), i.e. ^f T)S (t) = 
H s , r (0,t). In [PT3] the following crucial theorem is proved. 

Theorem 2.4. If k n > 0, then H TiS (w,t), H(w,t) and 

P w (t) = H(w,t)*{ty 1 
are polynomial functions on the variable t whose degrees are 

n-l 

deg H r>s (w,t) = w + y~]mm{rj,Si}, 

(12) . i=l 
deg H (w, t) = w + k\ — k n , 

deg P w (t) = w. 

It is important to point out that {P w } w >o is a sequence of matrix or- 
thogonal polynomials with respect to a matrix weight function W = W(t) 
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supported in the interval (0,1) and given in |PT3j . From (jlip it easily 
follows that {P w } w >o satisfies the following three term recursion relation 



(13) 



tPw{t) — A w P w -]_(£) + B w P w {t) + C w P w+ \(t). 



The above three term recursion relation which hold for all w > can be 
written in the following way 



(14) 



^0 




B 


Co 


Pi 






Bi 


P2 







A 2 


P 3 










Now we observe that the semi-infinite matrix M on the right hand side 
is a stochastic matrix, i.e. all the entries are nonnegative and the sum of 
the elements in any row is one. In fact, the elements in the r row of the w 
blocks are either zero or (A W ) T;S , (B w ) T>a , (C W ) T)S which are given in (fTUf) . 
Their sum is 

y~XA/;)r,s + (B w ) r , s + (C^) r>s = a£ +1 (m)&i (m + e n+ i) 



n+1 



+ a|(m)6?(m + e^) + ^ a] (m)&|(m + e,) 

i=2 j=i 

n n 

+ ^2 a 2 {m)b 2 n+l {m + ej ) + a 2 n+l (m) ^ b 2 (m + e n+1 ) 
+ ^ a|(m)6?(m + ej) + a\{m)b 2 n+l (m + e x ) 

2<i^j<n 

n 

+ a?(m)^&J(m + ei), 

i=2 

where we replaced m(w, r) by m. The right hand side can be rewritten to 
obtain 

n+1 

^2(A w ) r , s + (B w ) rja + (C w ) r>a = c£ +1 (m) ^ b 2 (m + e n+ i) 
sen j=i 

n n+1 n+1 

+ ^ a 2 (m) ^ b 2 (m + e,) + o?(m) ^ & n+ i(m + ei) 

j=2 i=l j=l 

n+1 n+1 

j=i i=i 
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Now by using ([6]) the assertion 



^ ^ (^w)r,s (Bw)r,s + (Cui)r,s — 1 



sen. 



follows, proving that the semi-infinite matrix M is stochastic. 



3. The substeps of the random walk 

In what follows we will construct a factorization of the stochastic matrix 
M appearing in (fl¥|) into the product of two stochastic matrices of the form 



(15) 



M 



Y X 

Y x X x 

Y 2 X 2 

o y 3 o 



5 o 

Ri Si 

o R 3 s 3 



While the random process given by the matrix M leaves invariant the set 
P introduced below, see (28), this is not true for its substeps going along 
with this factorization. This section deals with this complication in great 
detail. 

The multiplication formulas given in Theorem 12.21 restricted to g = a{6) 
give 



n+1 



(16) 



cos(W k (a(#)) = ^a2(m,k)$ s m+e - k (a^)), 

n+1 

co S (8)^ k (a(9)) = ^62(m,k)cl>r ej ' k (a(e)). 



We recall that we fixed k with k n > and we took m = m(w, r) as in 
([9|). Also making the change of variables t = cos(6>) we defined F rtS (w,t) = 

$s w,r ^' k (a(0)). Now we make the following important observation 



(17) m(w, r) ± ej = < 



m(w ± l,r) ±e n+ i if J = 1 , 

m(w, r ± ej_i) ± e n+1 ifj = 2,...,n, 
m(w,r) ± e n+ i ifj = n + l. 

Introduce the following scalar functions 

F+ s (w,t) = ^T {w ' r)+en+u \a(e)), 

and the matrix function 

F+(w,t) = (F+ s (w,t)) {rtS)enxn . 
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Then the first identity in (I16p becomes 

ViF riS (w,t) = al(m(w, r))F^ s (w + 1, t) 
(18) ^ , . , 2 
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+ ^ a ]+i( m ( w ^ r )) F ^+e 3 , s {w,t) +al +l (m(w,r))F^ s (w,t). 

3=1 

For each w>0we define the following matrix of type f2 x Q 



(19) 
where 



Xy; — ((^ui)r,s)) ^u) — ((^ / ii))r,s)) 

I af (m(u>, r)) if s = r, 



(^M>)r,s 



otherwise, 

«n+i( m O, r )) ifs = r, 
Oj +1 (m(w,r)) ifs = r + ej, 







otherwise. 



Now the set of scalar identities (|18|) with (r, s) E f2 x can be written as 
a matrix identity in the following more convenient way 

(20) V~tF(w,t) = X w F + (w + l,t) +Y w F+(w,t). 

For each w>0we define the following matrix of type f2 x Q 



(21) 
where 



Rw — \\Rw)r,s) j S w — ((5'io)r,s); 
J bl(m(w,r) + e n+ i) if s = r, 





(h2 



otherwise, 



b 2 n+l (m(w, r) + e n+ i) if s = r, 
b 2 j+1 (m(w,r) +e rt+ i) if s = r - e 







otherwise. 



If we multiply (|2U|) by \/i and use the second multiplication formula given 
in (I16|) we obtain 



(22) 



£F(w, t) =X w (i? w+ iF(u», t) + S w+1 F{w + 1, i)) 
+ y u ,(i?„,F(u; - 1, t) + S w F(w, t)) 
={X w R w+l + Y w S w )F(w, t) + X w S w+l F(w + 1, t) 
+ Y w R w F(w-l,t), 

since we claim that 

(23) VtF + (w, t) = R w F(w -l,t) + S w F(w, t). 
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Indeed we have 

VtF+ s (w, t) =Vt^T M+en+1 ' k (a(9)) 

n+l 

- m(w),r)+e n _|_i— ej,k , 



: £ b](m(w, r) + e n+1 )^' rJ+e " +1 - e - k (a(0)) 



=&5(m(u,,r) +e ri+1 )C (w - 1 ' r) ' k (a(9)) 



3=2 

+ b 2 n+1 ( m (w,r)+e n+1 )<l>? {w > r) ' k (a(e)), 

where we used (fT7|) . 
On the other hand 

{R w F{w - 1, i)) r>s = ^(i^) r , q F qiS (u; - 1, t) 

=bl(m(w, r) + e n+1 )F r , s (w - 1, t), 

and 

(S w F(w, t)) T)S = }XS w )r t(l Fq tS (w, t) 

=t% l+1 (m(w,r) +e n+1 )F riS (w,t) 

n-l 

+ ^ bj+i( m ( w > r ) + e„+i)F r _ ej>s (w, t). 

3=1 

Then fl23]) follows easily. 

Finally if we compare (|22p with (jlip in Theorem 12.31 we obtain 

which is equivalent to the factorization (|15p . 

We end by checking that both matrices in the right hand side of (|15p are 
stochastic: 

^ ^ (^w)r,s ~i~ ^ ^ (^Qo)r,s 

sen sen 

= a 2 n+1 {m(w,r)) + ^ af +1 (m(w, r)) + a?(m(w, r)) = 1, 
l<j<n-l 

^ ] (-^m)r,s "I - ^ ^ ('S l «))r,s 

sen sen 

= 6?(m(w,r) +e n+ i) + b 2 n+l {m{w, r) +e n+ i) 

+ 6 i+i( m ( u; ' r ) + e «+i) = !> 

l<i<n-l 

where we used that af (m, k) = 6?(m, k) = 1, see ©. 
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Now we want to consider the random walks associated to the probability 
matrices appearing in (fl~5j) . 



M 



B C 
A ± B x Ci 
A 2 B 2 C 2 



MiM 2 , 



(24) Mi 



X 
Yi Xi 

Y 2 X 2 



M 2 



S 
i?i S x 
i? 2 S 2 



Let F w and denote, respectively, the polynomial functions F u 
and = F+(t). Then ([23]) can be written as follows 



F w (t) 



(25) \/t 
Similarly (|20j) gives 

(26) Vi 



Ft 
Ft 



So 
Ri St 
R 2 S 2 



Fo 
F 2 



Fo 




Y 


Xo 







Fi 







Y\ 


x x 





F 2 









Y 2 


x 2 



Fo + 
Ft 

F+ 



We can now rewrite (I22p in matrix form, 



(27) 



*2 
Fo 



M 1 M 2 



Fo 




Fo 


Fi 


= M 


Fi 


F 2 


F 2 



The state space of the random walks W, W\, W 2 associated, respectively, 
to the stochastic matrices M, Mi, M 2 is the set N>o x f2, and W is equal to 
the composition W x o W 2 . 

We recall that the map (w, r) i-> m(u;, r) defined in ([9]) is an injection of 
N>o x into the k-spherical dual U(n + l)(k) of U(n + 1), and its image is 



(28) 



P = {m G U(n + l)(k) : s m = s k }, 



where s m = mH h m n+ i, s k = &i H h k n . 

Let us now consider the random walk W x associated to the stochastic 
matrix M x . Below we display the entries of Mi at the different sites of its 
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(w, r)-row, 



a 2 j+1 (m(w,r)) 



al(m(w,r)) 




if m(w, s)-site = m(u>,r), 
if m(u7, s)-site = m(w, r + e.,-), 
if m(w, s)-site = m(w + 1, r), 
in other sites. 



The appearance of the plus sign in the right hand side of (|26p makes it 
natural to consider instead the random walk obtained from W\ by 
applying a shift by e n+ i. Thus, if the system is at state m(w,r) at time t, 
then at time t + 1 it can move in the following ways 



with probability a^ l+1 (m(w, r)), 



r m(w,r) ->• m(w,r) + e n+1 , muup^a^j 

m(w, r) — > m(w, r) + ej+i, with probability a^ +1 (m(u>, r)), 

m(w, r) — > m(w, r) + ei, with probability a\ (m(w, r)), 

i m(u>,r) — > other states, with probability 0, 



because m(w,r + e^) + e n+ i = m(w,r) + e J+ i for 1 < j < n — 1, and 
m(w + 1, r) + e n+ i = m(?x;, r) + ei . This is in accordance with the following 
formula derived by looking at the ((u>, r), s)-entry of (|26p . 



n+l 



cos(#)$r Kr) ' k (a(#)) = 5^af(m(«;,r))$ 



m(ui,r)+e ;) - ,k 



(a(0)). 



Now it is worth to observe that does not leave invariant the subset 
P but extends to a random walk W\ in U(n + l)(k) defined by 

(29) Wx : m — y m + ej, with probability Oj(m, k). 

We proceed similarly with the random walk W2 associated to the stochas- 
tic matrix M2. Below we display the entries of M2 at the different sites of 
its (w, r)-row, 

if m(w,s)-site = m(w,r), 



'bl +1 (m(w,r) +e n+ i) 
b 2 j+1 (m(w,r) +e n+ i) 
&?(m(u;,r) + e n+ i) 




.1 1 ~ 



if m(w, s)-site = m(u>, r — e 
if m(w, s)-site = m(w — l,r), 
in other sites. 



The appearance of the plus sign in the left hand side of ([25]) makes it natural 
to consider instead the random walk W2 obtained from W2 by applying a 
shift by — e n+ i. Thus, if the system is at state m(w,r) at time t, then at 
time t + 1 it can move in the following ways 

'm(w,r) ->• m(tc,r) - e n+ i, with prob. b^ +1 (m(w, r) + e n+ i), 

w _ _ ^ m(w,r) -> m(w,r) - e j+1 , with prob. 6| +1 (m(w,r) +e n+ i), 

m(w, r) — >• m(w, r) — ei, with prob. 6f (m(u>, r) + e n+ i), 

k m(u;,r) — > other states, with prob. 0, 
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because m(w, r — ej) — e n+ i = m(w,r) — e J+ i for 1 < j < n — 1, and 
m(w — 1, r) — e n +i = m(w, r) — ei. This is in accordance with the following 
formula derived by looking at the ((w, r), s)-entry of (|25j) . 



n+l 



cos(6>)<£ 



m(to,r)+e n+ i,k 
s 



(a(0)) = ^>f(m(™,r)+e n+1 )<I> 

3=1 



m(iD,r)+e n+ i-ej,k, 



Then W 2 does not leave invariant the subset P but extends to a random 
walk W 2 in U(n + l)(k) defined by 

(30) W 2 : m — y m — Bj , with probability b'j (m + e n _|_i, k), 

for 1 < j < n + 1 . 

The transition matrices of W\ and PF 2 are, respectively, the following 
block bidiagonal matrices 



(31) Mi 



with 



Y X 
Yi X x 

y 2 x 2 



M 2 



5 
£1 5i 
R 2 S 2 



af(m) 


a n+l( m ) 

a| +1 (m) 




if n = m, 

otherwise, 

if n = m, 

if n = m + ej ; , 

otherwise, 



(-Rui)m.r 
(5»i>)m,r 



6?(m + e n+1 ) ifn = m, 

otherwise, 

6^ +1 (m + e n+ i) ifn = m, 
b 2 j+1 (m + e n+ i) ifn = r-ej, 




otherwise. 



where m, n G U(n + l)(k) are such that w = mi — k\ = n\ — k\, and 
1 < j < n - 1. 

Moreover, the stochastic matrix M corresponding to the composition 
W = W\ o W 2 is equal to M\M 2 , and it is given by 



M 



Co 
ii Bi C x 
A 2 B 2 C 2 
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with 



(A 



{Cu 



a 



-n+l( m ) & l( m + e 



i,n = < 



i,n = < 



a 



6 i+l 



(mjbf (m + e 



■n+l) 





«?( m )&n+l( m + e l)) 

a?(m)62 +1 (m + ei) 



i,n = < 



E 

l<j<n+l 



a 



,2 



a|(m)fe|(m + e,-)) 

( m ) & n+l( m + e j+l) 



if n = m 

if n = m + e,- 

otherwise 

if n = m 

if n = m — ej 

otherwise 

if n = m 



«n+l( m )& : 



, ,. ,^ +1 (m + e n+ i) 
;+i( m )^+i( m + e j 



m + 
m — 



e 3+l) 



if n 
if n 

if n = m + — 

otherwise 



where m, n 6 U(n + l)(k) are such that w = m\ — k\ = n\ — k\, and 
1 < i, j < n — 1. The coefficients a?(m), 6?(m) for 1 < i < n + 1 are those 
defined in ©. 

If we identify N>o x 17 with the subset P, defined in (|28p . by (u>,r) = 
m(zi;,r), then clearly W 
Therefore 



1¥|P, because M become a submatrix of M. 



W x oW 2 = W = W\ P = {Wi o Wb)|p. 

To conclude, the analysis of the random walk W associated to the stochas- 
tic matrix M is simplified by looking at the decomposition W = (W\oW2)\p 
instead of considering W = W\ o W 2 - 



4. An urn model for U(3) 

We now give a concrete probabilistic mechanism that goes along with the 
random walk W\ constructed in Section [3] by group theoretical means, see 
l\29\) . An entirely similar construction going with W<2 can be considered for 
the other substep of our process. 

This section is included for the benefit of the reader. It describes in 
detail, for the simple case of n = 2 going along with the pair (U(3),U(2)), 
a construction that will be given in general in Section [5j 

A configuration, or state of our system, is now a triple of integers m = 
(mi,m2,m3) subject to the constrains m\ > k\ > m,2 > k 2 > with 
two fixed integers k\ > k 2 , see (UJ. We describe a stochastic mechanism 
whereby one of the three values of the rrii is incresased by one with the 
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following probabilities, see ([5]) 




(m 1 -k 1 + l)(m 1 -k 2 + 2) 



ol (m,k) 



(mi - m 2 + l)(mi — m 3 + 2) ' 
(ki - m 2 )(m 2 — + 1) 




(mi - m 2 + l)(m 2 - m 3 + 1) ' 
(fei - m 3 + l)(k 2 - m 3 ) 



(mi - m 3 + 2)(m 2 - m 3 + 1) ' 



In the general scheme to be considered later this case corresponds to the 
value n = 2, and thus we start with two urns B\,B 2 . In urn Bj, j = 1,2, 
place rrij — kj + 1 balls of color c,- and fc,- — mi+i balls of color dj. These 
four colors are all different. Notice that we could have no balls of colors d\ 
or d 2 and that the total number of balls in urn Bj is rrij — mj+i + 1. 

It will be useful to consider the following ordered set of urns 



In view of the notation to be introduced in the general case we denote these 
urns as 



We will introduce later on an order among certain collections of urns that 
will yield, in this particular case, 



Now perform a total of three consecutive experiments. Each experiment 
consists of drawing one ball at random (i.e. with the uniform distribution) 
from an urn in the ordered set of urns above, record the outcome as a letter 
in a word, and continue to the next experiment making sure to return the 
ball that has been drawn to its original urn after this experiment has been 
performed. 

The first experiment consists of picking one ball from urn B\^\ = B\. This 
can give a ball of color c\ or d\. Record the outcome c\ or d\ as the first 
letter in a word of three letters, and return the ball to its original urn, B\ t \. 

The second experiment consists of picking one ball from urn B 2t2 = B 2 . 
This can result in a ball of color either c 2 or d 2 . Record the result as the 
second letter in a word that will have a total of three letters (the colors of 
the balls chosen in experiments 1,2,3), and return the ball to its original 
urn, £? 2 ,2- 

The last experiment consists of picking one ball from the union of the urns 
Bx i and -B 2>2 , i.e urn .Bi j2 . The color of the ball in question i.e. ci, d\, c 2 or 
d 2 is the last letter in our word. This last ball drawn from B\^ 2 = B\ U B 2 
is then returned to the urn B\ or B 2 where it came from. 

There is a total of sixteen (=2x2x4) possible words that can arise in 
this fashion from an alphabet of four letters. These words constitute the set 
of all possible outcomes of the experiment made up of these three succesive 
and properly ordered ones. 



Bi,B 2 , B\ U B 2 . 



B\ : i, B 2 2 , Bi j2 . 



#1,1 < B 2j2 < Bi i2 - 
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Since we return the chosen ball at the end of each one of these experiments 
to its original urn, we have that the state of the system has not yet changed. 
This is about to happen now. 

We need a rule to decide which of the three values mi,m 2 ,m3 will be 
increased by one unit as the result of our experiment. To this end we break 
up the set of sixteen words into three disjoint and exhaustive sets. These sets 
will be denoted by £1,3, £2,3 and £3,3, and the sample space £3 of cardinality 
16 is given by 

3 

Ss=U S jfl . 
j=i 

Each set £^3 consisting of words with three letters will be obtained by a 
"growth process" starting from the sets we would have if we had considered 
the previous case, namely n = 1, when we have only one box and we were 
dealing with U(2). In that case the sets are made up of words of one letter, 
either c\ or d\ . To make the connection with the general case we will denote 
these sets in the case of one urn by Sip and 6*2,2, and the sample space by 
£2 = £1,2 U £2,2. Explicitly Si i2 = {ci}\ S 2 , 2 = R}. 

Let us come back to the case n = 2. The class £13 is formed by including 
all three letter words that start as those of Si 2 and whose remaining two 
letters are such that the last one is not d 2 , i.e. either c±, di or c 2 . Thus 

Si,3 = {(ci,c 2 ,ci), (ci,c 2 ,di), (ci,c 2 ,c 2 ), (ci,d 2 ,ci), (ci,d 2 ,d 1 ), (ci,d 2 ,c 2 )}. 

The class £2,3 is formed by including all three letter words that start as 
those of £2,2 and whose remaining two letters are such that the first one is 
not d 2 . Explicitly £2,3 is 

S2,3 = {(di,c 2 ,ci), (di,c 2 ,di), (di,c 2 ,c 2 ), (di,c 2 ,d 2 )}. 

Notice that the meaning of the requirement "not d 2 is quite different 
when it applies to the second urn 7^2,2 as above, or to the third urn 73^2 as 
in the previous case. 

Finally £3,3 is obtained by taking the union of all three letter words that 
start as in Si^ and have d 2 as their last letter, together with all words that 
start as in £2,2 an d have d 2 as the second letter. Notice that £3,3 is obtained 
by going over all the classes already built, £1,3 and £2,3, and replacing the 
condition not d 2 by d 2 . The class £3,3 is thus made up of two sets of words, 
namely 

£3,3 ={(ci,c 2 ,d 2 ), (ci,d 2 ,d 2 )} 

U {(di,d 2 ,ci), (di,d 2 ,di), {di,d 2 ,c 2 ), (di,d 2 ,d 2 )}. 

It takes almost no effort to see that all these 6 + 4 + 6 = 16 words have 
been classified into three disjoint and exhaustive classes. 
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Now we compute the total probability of getting a result that belongs to 
each class. For the first class Si 3 we have, 



(m 1 -fci + l)(mi-A& + 2) 



o?(m,k). 



(mi - m 2 + l)(mi - m 3 + 2) 
For the second class 5 2i 3 we have that the probability is 

(ki -m 2 )(m 2 - k 2 + 1) 2 , i\ 

7 — — = a 2 (m, k). 

(mi - m 2 + l)(m 2 - m 3 + 1) 

Finally the total probability of the third class £3,3 is, 

(mi - ki + 1)(Aj 2 - m 3 ) (&i - m 2 )(A; 2 - m 3 ) 



+ 



(mi - m 2 + l)(mi - m 3 + 2) (mi - m 2 + l)(m 2 - m 3 + 1) 
_ (fc 2 -m 2 )(A;i -m 3 + l) _ 2( , 
(mi - m 3 + 2)(m 2 - m 3 + 1) 

We are ready to give a rule for changing the state of the system in one 
unit of time. A result belonging to the subset Sj t s, j = 1,2,3, will lead to 
a transition to a new state m + ej, where rrij is increased by one. In terms 
of balls this will be achieved by removing from each urn containing a ball 
of color dj-i one of these balls, and adding to each urn containing a ball of 
color Cj one ball of this color from the bath. When j = 1 we do no removal. 

5. AN URN MODEL FOR EVERY U(n + 1) 

In this section we describe a random mechanism that gives rise to a 
Markov chain whose one-step transition matrix is 

Y X 
Yy X x 

Y 2 X 2 

o y 3 * 3 o 

appearing in the factorization (|15[) and where the matrices Xj, Y{ are defined 
in (HID. 

A configuration is a set of n + 1 values of the integers rrii, 1 < i < n + 1, 
subject to the constrains mi > &i > m 2 > • • • > m n > k n > m ra +i where 
the integers fe, remain unchanged throughout time. We will construct a 
stochastic process whereby in one unit of time one of the rrij is increased by 
one with probability given by 



(32) 



(m,k) 



Yiiftimi - rrij - i + j) 

Consider n urns Bi, . . . , B n . In urn Bj place rrij — kj + 1 balls of color 
Cj and kj — rrij + i balls of color dj. We assume that the colors Cj,dj are all 
different. Notice that in urn Bj may be no ball of color dj, and that the 
total number of balls in Bj is rrij — rrij+\ + 1. 
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Consider the following ordered set of urns 



B\, B2, B1UB2, -B3, B2UB3, B1UB2UB3, . . . , B n , B n _iUB n , . . . , BiU- ■ -UB,, 



The union of urns is an urn whose content is the union of the set of balls 
in each urn in the union. Observe that the total number of urns under 
consideration is n(n + l)/2. Let 



is ordered lexicographically according to: (k,j) < (r, s) if j < s or if j = s 
and r < k. 

We will perform a total of n(n + l)/2 consecutive experiments. Each 
experiment consists of drawing one ball at random (i.e. with the uniform 
distribution) from each urn in the ordered set of urns, record the outcome 
as a letter in a word, and continue to the next experiment making sure to 
return the ball to the original urn after this experiment has been performed. 
One should think of a complete experiment as consisting of these n(n + 1) /2 
individual experiments. The transition from the present state of the system 
to the next one takes place after the complete experiment is carried out. 

The first experiment consists of picking one ball from urn B\ t i, this can 
give a ball of color c± or d\. The result is recorded and the ball is put back 
in urn B\ \. The second experiment consists of picking one ball from urn 
#2,2 , this can result in either a ball of color C2 or 6,2- Record the result as 
the second letter in a word that will have a total of n(n+ l)/2 letters. Put 
the ball back in urn -62,2- Keep on going by taking successively at random a 
ball from an urn B k j and adding the letter corresponding to its color to the 
right of the word obtained in the previous step. The process finishes once a 
ball of the last urn B\^ n is picked and a final word of n(n + l)/2 letters is 
obtained. 

The alphabet is the set {cj, dj : 1 < j < n} of 2n letters. Then the sample 
space S n+ i consists of all words w of n(n + l)/2 letters that can be written 
with such an alphabet with the restriction that the letters allowed in the 
place (k,j) correspond to the color of any ball in urn B^ j. The cardinality 
of the sample space is 



Now by induction on n > 1 we define a partition of S n+ \ into n+1 disjoint 



B ktj = B k UB k+1 U---UBj, l<k<j. 
Clearly Bjj = Bj, and the set of all urns 

{B kJ : 1 < k < j < n} 




n 2(j~k+i). 



l<k<j<n 



subsets 



n+1 




3=1 
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For the benefit of the reader the construction will be spelled out in detail 
for small values of n after we describe it in the general case and prove 
Proposition 5.2. 

We start with S 2 = Sip U 6*2,2 where 

5*1,2 = {d 1 }, 5*2,2 = {di}, 4i=ci. 

Then 

15*1,21 = \^2,2\ = 1, 1 5*2 1 = 2- 

We make the following convention: the symbol <ij in the (fc,j)-place of 
a word stands for any color of a ball in urn B^j different from dj, and the 
letter x in the (k, j)-place of a word stands for any possible color of a ball 
in urn B^j. 

If n > 2 we set 

5l, n +l = {Wl,n+1 = Wi,nX ■ ■ ■ xd n G S n+1 \ W^ n G S*i,„}. 

Observe that the number of letters in the word iOi, n +i to the right of the 
word wi^ n is n. Similarly we define 

52,n+l = {w 2 , n +l = W 2 , n X ■ ■ ■ xd n X G S n+ i 1 W 2 , n G 5 2 ,„}. 

More generally for 1 < j < n we let 

Sj,n+1 = {Wj,n+1 = Wj,nX ■ ■ ■ X d n X ■ ■ ■ X G S n+ l : U>j,„ G S*j, n } 

where the number of x's to the right of ^ n is j — 1. 
The definition of 5 ra+ i, n+ i is more interesting, namely 

5 n +l, n +l ={Wn+l,n+l = Wi,nX ■ ■ ■ xd n G 5 n+ i : Wi jTl G 5i, n } 

U {w n+ i, n+ i = w 2 ,n^ • • • xd n x G S n+ i : W 2 ,n G S*2, n} 

U • • • U {u; n+ i )n+ i = w n;n d n x ■ ■ ■ x G S* n+ i : iw n , n G 5 n , n }. 

Proposition 5.1. Let n > 2. T/ien /or 1 < j < n we have 

\S hn+ i\ = \S j>n \{2{n-j) + l) 11 2(n-fc + l), 

l<fc<n, fc^j 

|5 n+ i, n+ i| = ^ |S*j, n | ]| 2(n-fc + l). 

l<?<ra l<fc<n, fc^j 

Proposition 5.2. {5j, n+ i : 1 < j < n + 1} is a partition of the sample 
space S n+ i. 

Proof. The proof is by induction on n > 1. For n = 1 we have 

5*2 = {d^ di}, Sip = {di}, 5 2 , 2 = {di}. 

Thus the statement is true for n = 1. Now assume that S n = Uj=i 5j, n 
is a partition of S* n for n > 1. If u> G 5 n +i, then w; = Wj >n x ■ ■ ■ x where 
Wj, n G S*j, n for a unique j. The x in the j-place of the last n letters is either 
d n or of the form fi n . In the first case w G S* n+ i, ra+ i and in the second 
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case w G Sj <n+ i. Thus S n+ \ = Uj=i Sj,n+i- At the same time we saw that 
w G £j,n+i for a unique 1 < j < n + 1. This completes the proof. □ 

The construction above is now made explicit for small values of n. 

1) n = 2. 

<Si,3 = {{iiXfi 2 }, S 2 ,3 = {d 1 d 2 x}, S 3t3 = {d 1 xd 2 } U {did 2 x}, 
|Si, 3 |=6, 152,31=4, 1^3,3] = 6, |S 3 | = 16. 

2) n = 3. 

<Si,4 = {fax fi 2 xx ■> ^2,4 = {did 2 xxd 3 x}, 
53,4 = {d 1 xd 2 d 3 xx} U {diG^x^gXx}, 
54,4 ={^1X^22^3} U {did 2 xxd 3 x} 

U {^iXf^o^xx} U {did 2 xd 3 xx}, 
\Si,i\ = 240, |S 2>4 | = 144, |5 3 , 4 | = 144, |5 4 , 4 | = 240, |5 4 | = 768. 

3) n = 4. 

5l,5 = {^iX^Xxd^XXxd^, 52,5 = {^1 ^ 2 a;a: '^3 xxx ^4 x }' 

^3,5 = { j^j xc^ fi-^xxx d^xx} U {d\d 2 x d 3 xxx d 4 xx} , 
$4,5 ={d 1 xd 2 xxdsd 4: xxx} U {did 2 xxd^xd^xxx} 

U {dixd 2 d$xx d 4 xxx} U {c/i d 2 xd-$xx d^xxx}, 
^5,5 ={$iX d 2 xxd 3 xxxdi} U \d\d 2 xxd 3 xxxdt±x\ 
U {d 1 xd 2 d 3 xxxd4 : xx} U {did 2 x ji^xxxd^xx} 
U {dix ^xxd^d^xxx] U {eZi d 2 xxd 3 xd^xxx} 
U {dxxdqdzxxd^xxx} U {did^xd-ixxd^xxx} , 
1^1,5 1 = 80640, | £2,5 1 = 46080, |5 3>5 | = 41472, 
|5 4 , 5 | = 46080, |5 5 , 5 | = 80640, |5 5 | = 294912. 



Theorem 5.3. T/ie probability to obtain a word w G Sj,n+i *s Oj( m )k) /or 
aii 1 < i < n + 1. 

Proof. Given (m, k) let m' = (mi, . . . , m n ) and k' = (Ari, . . . , k n -i). Then 
from (j32|) we get 

a (m,k) = a (m ,k) ■ — -, 

J J rrij - m n+ i + n - j + 1 

for all 1 < j < n. This result allows us to prove the theorem by induction 
on n > 1. When n = 1 we have only one urn B\ with mi — k\ + 1 balls of 
color ci and fci — m 2 balls of color cii . Thus the probability to obtain a word 
in Si 2 is 



mi — ki + 1 -, 



mi — m 2 + 1 



Oi(m,k), 
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where m = (mi, 7712) and k = (ki). Similarly the probability to obtain a 
word in £2,2 is 

h-m 2 2 

= a 2 (m, k). 

mi — m,2 + 1 

Thus the theorem holds for n = 1. Now assume that the theorem is true for 
n > 1. If 1 < j < n we have 

-Sj>+1 = {^j,n+i = ^j,n^ • • • x^ n x • • • x G 5 n+ i : itfj in G S^n} 

where the number of x's to the right of $ n is j — 1. Thus the probability 
to obtain a word w G S^n+i is equal to a| (m^k') times the probability to 
obtain the symbol $ n from the urn Bj )U . Now we recall the composition of 
urn Bj )Tl . By definition 

Bj j7l = Bj U Bj + i U • • • U B n , 

the total number of balls \Bj n \ = rrij — m n+ i + n — j + 1 and the number 
of balls of color d n is fc n — m n +i. Therefore the probability to obtain the 
symbol $ n from urn Bj )Tl is 

rrij — k n + n — j + 1 



to, - m n+ i + n - j + 1 
Hence the probability to obtain a word w G •Sj.n+i is 

a Am ,k) ^ ■ — - = od(m,k), 

J rrij - m n+ i + n - j + 1 J 

which establishes the theorem for all 1 < j < n. Since ^i<j< n +i a j ( m ' ^) = 
1 (see ([6])) and ^n+i = Ui<j< n +i Sj,n+1 is a partition of S n+ \ it follows that 
the statement of the theorem is also true for j = n + 1. □ 

Since we return the chosen ball at the end of each individual experiment 
to its original urn, we have that the state of the system has not yet changed. 
This is about to happen now. 

The outcome of a complete experiment produces a word that belongs to 
one of the subsets S^n+i in the partition of the sample space S n+ i. De- 
pending on which subset turns up we take a different action, thus obtaining 
a random walk in the space of configurations m = (mi, . . . , m n +i) which 
satisfy the constraints mi > k\ > ■ ■ ■ > m n > k n > m n +i imposed by the 
fixed n-tuple k = (ki, . . . , k n ). This simple process will give for each config- 
uration m a total of at most n + 1 possible nearest neighbours to which we 
can jump in one transition. 

A result belonging to the subset Sj )n +i, j = 1, . . . , n + 1, will lead to a 
transition to a new state m + ej , where rrij is increased by one. In terms 
of balls this will be achieved by removing from each urn containing a ball 
of color dj-i one of these balls, and adding to each urn containing a ball of 
color Cj one ball of this color from the bath. 
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Notice that all these transitions keep the values of k\ , . . . , k n unchanged 
and any transition that would violate the constrains does not occur because 
the corresponding probability a"-(m, k) vanishes. 



To each configuration m\ > k\ > m 2 > • • • > m n > k n > m n +i > we 
associate its Young diagram which has mi boxes in the first row, k\ boxes 
in the second row, and so on down to the last row which has m n+ \ boxes. 



We will construct a stochastic process whereby in one unit of time one of 
the mi is increased by one with probability a?(m, k) see ([5]). As in Section 
[5] this will require running some auxiliary experiments. 

We start with the case n = 1. We perform the following experiment to 
decide if we will increase m\ or m^- we choose to insert a box among one of 
the mi — k\ last boxes of the first row or to delete a box from the k\ — m2 last 
boxes of the second row. An insertion can occur either to the left or to the 
right of one of the m\ — k\ last boxes. We observe that there are m\ — k\ + 1 
possibilities of an insertion and k\ — m 2 possibilities of a deletion. All these 
are assigned the same probability. 

As an output of the experiment we get either a diagram with m\ + 1 boxes 
in the first row, or a diagram with k\ — 1 boxes in the second row. Here we 
are implicitly assuming that k\ > mi- If k\ were equal to m 2 we would get 
no Young diagram. Thus the sample space S of our auxiliary experiment 
consists of two (or one) Young diagrams which are obtained from the original 
one by adding one box to its first row or deleting one from its second row. 
Let Si be the subset of S consisting of the diagram with one more box in 
the first row, and let S 2 be the subset of S consisting of the diagram with 
one less box in the second row (or the empty set). Then the probability to 
obtain a diagram in S\ after the experiment is performed is 



6. A Young diagram model for U(3) 



D = 



Figure 3. m = (8, 



5,1), k 



(6,3). 




Similarly the probability to obtain a diagram in S2 is 

k\-m 2 , , , 2 
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as we wished. In the first case we go from the state (m, k) to (m + ei, k), 
and in the second case we go from the state (m, k) to (m + e2, k). 

Now let us assume that n = 2. In this case we will perform three con- 
secutive auxiliary experiments. The first experiment consists of inserting a 
box among one of the m\ — k\ last boxes of the first row or of deleting a 
box from the k\ — 777 2 last boxes of the second row. The second experiment 
consists of inserting a box among one of the 777,2 — &2 last boxes of the third 
row or of deleting a box from the fc 2 — 7773 last boxes of the fourth row. 
Finally the third experiment consists of inserting or deleting a box in one 
of the first four rows of the diagram as we did in the previous experiments; 
odd rows go along with insertion and even rows with deletion. If k\ > 777 2 
and &2 > 7773 the complete experiment gives rise to a triple (Di, D2, -D3) of 
Young diagrams: D\ is obtained from the original one by adding one box 
to its first row or by deleting one box from the second row, D2 is obtained 
from the original one by adding one box to its third row or by deleting one 
box from the fourth row, and D3 is obtained by adding one box to the first 
or to the third rows of the original diagram or by deleting one box from the 
second or the fourth rows. 

In what follows we use the following notation: D denotes the Young 
diagram corresponding to the original configuration (m, k) and D' = D±ej 
denotes, respectively, the diagram obtained from D by adding or deleting one 
box to the j-row of D, j = 1, 2, 3, 4. Observe that the sample space consists 
of all triples of Young diagrams (Di, D2, -D3) with D\ = D + ei, D — e2, 
D 2 = D + e 3 , D - e 4 , and D 3 = D + ei, D — e 2 , D + e 3 , D - e 4 . 



D + ei = 















+ 















































Figure 4. m = (9, 5, 1), k = (6, 3). 



D -e 2 = 



Figure 5. m = (8, 5, 1), k = (5, 3). 
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Thus our sample space S3 has generically 2 x 2 x 4 = 16 elements. The 
cardinality of £3 can be smaller, for example if k± = m 2 and k 2 7^ m 3 , then 
|S 3 | = 6. 

Let us partition the sample space S 3 into the following three classes. 
(33) 

5i, 3 ={{D X ,D 2 , D 3 ) :D 1 = D + e 1 ;D 2 = D + e 3 ,D-e A ; 

D 3 = D + e 1 ,D + e 3 ,D-e 2 }, 
S 2 ,3 ={(D 1 ,D 2 ,D 3 ) : D 1 = D-e 2 ;D 2 = D + e 3 ; 

D 3 = D + ei, D + e 3 , D - e 2 , D - e 4 }, 
S 3)3 ={(I>i, D 2 , D 3 ) : D x = D + ei; D 2 = D + e 3 , D - e 4 ; D 3 = D - e 4 } 
U {(L>i, D 2 , D 3 ) : D 1 = D - e 2 ; D 2 = D - e 4 ; 

D 3 = D + ei, D - e 2 , D + e 3 , Z? - e 4 }. 

We have |Si j3 | = 6, |S 2;3 | = 4 and |S 3)3 | =2 + 4 = 6. By simple inspection 
we see that S3 is the disjoint union of Si ;3 , S 2)3 and S 3;3 . 

Then the probability to obtain a diagram in Si i3 after a complete exper- 
iment is performed is 

(mi - fci + 1) (mi - fc 2 + 2) = fl2 
(mi-m 2 + l)(mi-m 3 + 2) U ' 
Similarly the probability to obtain a diagram in S'2,3 is 

- (fcl " m2) - ^ ~ k2 + '\ = a\ (m, k) . 
(mi —7712 + 1) (m 2 - m 3 + 1) 

Finally the probability to obtain a diagram in £3,3 is 

(mi - ki + 1) (mi —7712 + I) (fci - 7712 ) (mi - m 2 + 1) 



(fc 2 - m 3 ) (mi - m 3 + 2) (fc 2 - m 3 ) (m 2 - m 3 + 1) 

o§(m,k), 



(fe 2 - m 2 )(k 1 - m 3 + 1) 



(mi - m 3 + 2)(m 2 - m 3 + 1) 
as desired. 

If k\ = m 2 and k 2 7^ m-3 then |Si )3 | = 4, S 2)3 = and |S 3;3 | = 2. The 
probability to obtain a diagram in Si )3 is 

mi - k 2 + 2 2 
= af(m,k). 

mi-m 3 + 2 iV ; 

The probability to obtain a diagram in S 2i3 is = a 2 (m, k), and the proba- 
bility to obtain a diagram in S 3)3 is 

k 2 -m 3 2 

— = a 3 (m,k), 

mi — m 3 + 2 

as expected. 

Now the state of our random walk is modified in one unit of time as 
follows: if the outcome of the complete experiment above belongs to S Ji3 , 
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then we go from (m, k) to (m + e», k), j = 1,2, 3. In terms of diagrams we 
move from D to D + e2j-i, j = 1, 2, 3. 

7. A Young diagram model for every U(n + 1) 

Given a Young diagram D corresponding to the original configuration 
(m, k), D' = D db ej denotes, respectively, the diagram obtained from D 
by adding or deleting one box to the j-row of D, j = 1, . . . , 2n + 1. The 
stochastic process we are going to construct will have a transition mechanism 
determined by first performing a sequence of auxiliary experiments E^ j to 
be described now. We start by considering the following set of consecutive 
pairs of rows of the diagram D, 

{(l,2),(3,4),...,(2n-l,2n)}. 

The experiment Ekj, 1 < k < j < n, consists of inserting at random a 
box in an odd row i among the last mi — ki last boxes of such a row, or 
deleting at random a box in an even row i from the last ki — m^i last boxes 
of such a row. The row i is also chosen at random in the set of consecutive 
rows 

{2k- 1,2k,..., 2j}. 

The sequence of experiments is obtained by ordering them by the lexico- 
graphic order E^j < E r>s if j < s or j = s and r < k. Thus our sequence is 
the following one 

E2,2, E\p, -^3,3) -^2,3, -Ed,3j • • • > E n ^ n , E n _i :n , . . . , E\^ n . 

The symbol D ± ^ in the place corresponding to the experiment E^j of 
an ra(n + l)/2-tuple of diagrams, will stand for any possible outcome of E\. j 
except the diagram D ± e^, respectively. While an X in such a place stands 
for any outcome of E^j. For example in the case n = 2 considered before, 
see (I33p . we can write 

5i, 3 ={(D- &,X,D- 
S 2 , 3 ={(D-e 2 ,D- &,X)}, 

S 3 , 3 ={(D - X,D — e 4 }) U {(D -e 2 ,D- e 4 , X)}. 

Now we have a convenient notation to define inductively, for n > 2, a 
growth process similar to the one introduced in Section [5j to break up the 
outcomes of the sample space S n+ i into sets S^n+i (j = 1, . . . , n+1) starting 
from the partition of S n into sets Sj >n (j = 1 , . . . , n) . Let D^ n denote any 
n-tuple in the set <Sj jn , then we set 

S\,n+\ — \P\,n+\ = (Di :n ,X, • • • , X, D — $ 2n ) G S n+ i : Di jU G Si iV }. 

Observe that the number of diagrams in the (n + l)(n + 2)/2-tuple D% !n+ i 
to the right of the n(n + l)/2-tuple D\ >n is n. Similarly we define 

5'2,n+l = {E>2,n+l = (-t>2,n, X, - ■ ■ , X, D — $ 2n , X) G -S n+ i : D 2 ^ n G S2,n}- 
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More generally for 1 < j < n we let 

Sj,n+1 

= {Dj >n +i = {Dj, n , X, - ■ ■ ,X,D — ^2mX, • • • , X) G S n+ \ : Dj :Tl G Sj >n } 

where the number of X's to the right of D — ^ 2 n is J — 1- 

The definition of S n+ i :n+ i is (as in Section [5]) more interesting, namely 

S n +l,n+l = { Ai+l,n+l = (Di : n, X, • • • ,X,D — e2 n ) S S n+ i : Di :n G Si iU } 

U {-D n +l,n+l = {D2,n, X, - ■ ■ ,X,D — e 2n , X) G SVj+1 : Z?2,n G S^n} 

U • • • U {D n+ i^ n+ i = {D n , n , D — e 2n , X, • • • , X) G 5 n+ i : -D nn G S n ^ n }. 

Now by induction on n > 2 it is easy to prove that {Sj jn +i : 1 < j < n+1} 
is a partition of S" n+ i. Also by induction on n > 2 it is possible, as we did 
to established Theorem 15.31 to prove the following main result. 

Theorem 7.1. The probability to obtain an n(n + 1) /2-tuple of diagrams 
Dj^n+i £ Sj,n+i is a?j (m, k) (see (j32j) ) for all 1 < j < n + 1. 

The outcome of a complete experiment produces an n(n + l)/2-tuple 
of Young diagrams that belongs to one of the partition subsets <Sj>+i of 
the sample space S n+ \. Depending on which subset turns up we take a 
different action, thus obtaining a random walk in the space of configurations 
m = (mi, . . . , m n +i) which satisfy the constraints 

mi > k% > ■ ■ ■ > m n > k n > ra n+ \ > 0, 

imposed by the fixed n-tuple k = (ki, . . . , k n ). This simple process will give 
for each configuration m a total of at most n + 1 possible nearest neighbours 
to which we can jump in one transition. 

A result belonging to the subset Sj }Tl +i, j = 1, ... ,n + 1, will lead to a 
transition to a new state m + e^, where rrij is increased by one. 

Notice that all these transitions keep the values of k\ , . . . , k n unchanged 
and any transition that would violate the constrains does not occur because 
the corresponding probability a^(m, k) vanishes. 
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